How to mode-lock an atom laser 
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Abstract 



We consider a possible technique for mode-locking an atom laser, based on the gen- 
eration of a dark soliton in a ring-shaped Bose-Einstein condensate, with repulsive 
atomic interactions. The soliton is a kink, with angular momentum per particle equal 
to h/2. It emerges naturally, when the condensate is stirred at the soliton velocity, 
. and subsequently outcoupled with a periodic Raman pulse-train. The result is a repli- 

cating coherent field inside the atom laser. We give a numerical demonstration of the 
generation and stabilization of the soliton. 

> 

The discovery of Bose-Einstein condensation (BEC) in ultra-cold alkali-metal vapors |lj in 
a magnetic trap, at temperatures of lO -6 ^ or lower, has raised the possibility of a coherent 
atom laser. Recent experimental developments show that this is indeed practical. All 
atom lasers to date, however, produce output pulses that are not phase- coherent. It would 
be desirable to have a mode-locked laser, in which the pulses are in phase with each other, for 
this will enable a wide range of interference experiments and phase-sensitive measurements. 
In addition, mode-locked lasers can have an enhanced intensity stability, relative to their 
non-mode-locked cousins. 

One way to make a mode-locked atom laser is to create a periodic field circulating around 
a ring-shaped condensate, and outcouple it with a synchronized period. To do this, we must 
choose a periodic field that can be easily created, and that has sufficient stability to enable 
a steady-state outcoupling process. In this regard, we suggest a dark soliton in a condensate 
of atoms with repulsive interactions. This is a kink configuration, which stands apart from a 
continuum of possible excitations, owing to two distinctive features: (a) It has a characteristic 
propagating velocity v ; (b) The angular momentum per particle, normal to the plane of the 
ring, is h/2. 
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The dark soliton can be created by stirring the condensate at the characteristic velocity 
Vo] but it has to be cleaned up, because the stirring also creates other excitations such as 
phonons. The cleansing can be achieve by applying a stroboscopic loss mechanism, which also 
serves as outcoupler for the atom laser. The soliton persists because of a dynamical stability 
due to the stroboscopic environment, which is particularly strong for the half- integer angular 
momentum dark kinks. In the following, we first describe the soliton as a solution of the 
nonlinear Schrodinger equation (NLSE) and then demonstrate mode-locking via numerical 
simulations. 

Consider a condensate contained in a ring of radius R. Let the cross-sectional radius be 
r , and let the cross-sectional area be denoted by A = nr\. We choose R 3> r , so that the 
transverse excitations are far more energetic than those along the ring. For the low excitation 
modes, therefore, we may regard the condensate as a one-dimensional system, and denote 
by 9 the angle around the ring. The condensate wave function ty(9,t) salsifies the NLSE 

h 2 d 2 ^> TTT ^h 2 a^ aT 

%h — = hCAtH h&r^ (1) 

dt 2mR 2 dd 2 m 1 1 K ' 

where m is the atomic mass, a the s-wave scattering length, and U is an external potential 
due to trap non-uniformity. The total number of atoms N enters through the normalization 
ARJq d6\^\ 2 = N, which is a constant of the motion, if U is independent of the time. 
Continuity requires the boundary condition ty(6 + 2n,t) = ty(6,t). 

The NLSE on a line is well-known in nonlinear optics, where it describes the envelope of an 
electromagnetic wave propagating along an optical fiber ||. Let us recount what is generally 
known in the linear case. For attractive interactions with a < 0, the nonlinear term in the 
equation presents an attractive potential proportional to |\l/| 2 . In three spatial dimensions 
in free space, this would lead to "self-focusing" - - the development of spots of infinite 
intensity in finite time [|J. In one dimension, however, the kinetic energy counterbalances 
the attraction, leading to the formation of a stable bright soliton. The situation is the same 
when we join the ends of the line to form a ring. 

For repulsive interactions with a > 0, we can make a dark soliton in a linear condensate 
by requiring \1/ to approach ±1 at opposite ends of the line. By continuity, the configuration 
must have a kink, i.e. a zero of the wave function. The slope at the kink determines its 
propagation velocity vq, which will be finite and nonzero. Kink solitons in a cigar-shaped 
BEC were recently observed ||. Gray solitons also exist, in which the wave function has a 
minimum, but never vanishes. These have also been analyzed theoretically ||, and created 
experimentally J7|. They propagate with a speed proportional to the wave intensity at the 
minimum, and hence come to rest in the dark-soliton limit. 

These belong to a different class from the kink soliton we are considering on a ring, where 
periodicity demands that the ends match. Therefore, the phase of the wave function must 
change by nix upon one complete revolution, where n is an odd integer. This makes the 
angular momentum per particle normal to the ring nfi/2, for the same mathematical reason 
that an electron has spin 1/2. The lowest-energy dark soliton has n = 1. 
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It is convenient to use dimensionless variables to reduce the NLSE to the form 

dip d 2r i 
~8r ~ ~~89 



M = _tH +vi) + a ^ (2) 



where r = (h/2mR 2 )t, ip = ^^/2tcAR/N, V = {2mR 2 /h 2 )U, and a = ANRa/A. The 
boundary condition is ip(9+2-K, r) = ip(9, r), and the normalization condition is J Q 27r d9ip*ip = 
2tt. The energy per particle e, and angular momentum per particle i normal to the ring are 
given by 



d9 



dip* dip a , . 2 - 



Am J V 89 r 89 

We consider only special solutions relevant to the application at hand, and discuss more 
general solutions in a separate paper ||. 

For configurations that move around the ring at speed vq, we denote by 9 = 9 — vqt the 
angle in the co- moving frame, and put 

^(0, r) = f(9) exp (i£9 - iur) (4) 

where / is real, lu is a constant to be determined by solving the equation, and I is an integer 
or half-integer. The boundary condition requires f(9 + 2n) = (—l) 2i f(9). Substituting (ffl) 
into (^), we find that the angular momentum per particle is £: 

£ = 0, |, 1, §,--.. (5) 
Substituting / into (0), we find the propagation velocity 

v = 2£. (6) 

The equation for / is 

where j3 = u + I 2 . Mathematically, this equation describes the Newtonian motion of a 
particle of unit mass, in a potential $(/) = — f f 4 + § f 2 - From the conserved "energy" 
E = \f' 2 + $(/)> we can fi n< i / i n terms of elliptic functions. Gray solitons correspond to 
complex /, and are not considered here. 

From now on we assume a > 0. In the absence of external potential, a uniform state 
must have ip*ip = 1 by normalization. These are the vortex states: 

VWtex(0 , t) = exp [i£9 - i(a + £)t] 

evortex = 2+ £2 (^ = 0,1,2,3,...) . (8) 



4 



P.D. Drummond, A. Eleftheriou, K. Huang, and K.V. Kheruntsyan 



In contrast, the solutions with half-integer £ must vanish somewhere, and thus correspond to 
dark solitons. The width of the soliton should be proportional to iV~ 1//2 , for the correlation 
length of the condensate is (87rap) -1 / 2 , where p is the particle density. In the large N limit, 
we obtain the simple form: 



/soliton (0) ~ t&s&(J^(Q - it] 



where a oc # > 1. In this limit we have f3 ~ a + (27r) _1 v^2a. These limiting solitons all 
have the same shape, but propagate at different velocities 2£. Derivation of these and more 
general results will be given in a separate paper [|J. 

In experimental situations, the condensate is in equilibrium with a thermal cloud of 
uncondensed atoms. The equilibrium number of condensate atoms, relative to that in the 
thermal cloud, is determined by the temperature, and it determines the width of the soliton. 
During each soliton round-trip, a loss of particles from the condensate due to output coupling 
will be countered by a gain from the thermal cloud. The mean atom number in the entire 
BEC is self-stabilizing because the output-coupling efficiency is a nonlinear function of the 
dark-soliton width, which in turn is inverse to the total atom number. Thus, a fluctuation 
from equilibrium that increases the atom number, will also increase the efficiency of the 
output coupler, which therefore acts as a nonlinear absorber. Even with linear gain and loss 
mechanisms, we can create a self-maintained soliton and a steady stream of coherent output 
pulses. The feasibility of the scheme can be demonstrated by numerical simulations, as we 
now describe. 

First, we create a dark soliton with i = 1/2 from a uniform static condensate, by mo- 
mentary stirring it at the soliton velocity Vq = 1. In practice, this can be produced by a blue 
detuned laser beam. In the simulation, we introduce a repulsive potential V{9, r) to create 
a moving "hole" in the condensate at 9 = r — it. The time origin is displaced by it, to wait 
for the hole to fully form. Next, we have to clean up the configuration by adding gain and 
loss mechanisms, for the stirring creates other excitations in addition to the dark soliton. 
The entire procedure is contained in the generalized equation of motion 

l % = -f^r + «H 2 ^ + v (°> T W + i l9- ^ r )] ^ ( 10 ) 

where the stirring potential is chosen as: 

V(^) = ^exp[-(^±^) 2 -(I^) 4 ]. (u) 

In the gain-loss mechanism, g is a constant gain rate, representing stimulated emission from 
non-condensed atoms that are continuously loaded into the trap. The loss function makes 
localized periodic hits at the center of the soliton: 
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Figure 1: Modulus of the condensate wave function as function of angle around the ring, and 
dimensionless time. Stirring an initially uniform condensate at the soliton velocity creates a 
moving hole, which develops into a dark soliton, after unwanted excitations are cleansed by 
the loss mechanism. 



7(0, r) = 7oexp 



e \ 2 (t-t x 



- 



Oq) \ U 



(12) 



where T\ = (r mod 2n) — ir . This simulates a stroboscopic output coupler, realizable, for 
example, through a local Raman-tuned transition to a non-trapped state of the atoms. The 
strobe acts for the first time at r = 2n, allowing time for the hole to form. 

In the numerical analysis, we set a = 7.3634, chosen to correspond to /3 — 8 in (0). For 
the stirring potential we use V Q = 1.673a , cti = 1.575^27/3, and a 2 = tt/2. The gain and 
the loss parameters are set to: g = 0.01, 70 = 0.9, oq = 1.05y2//3, cr r = 1.05-^/2//?. 

Following are the results of numerical calculations, with the initial condition ip = 1. 
Fig. |l] shows j^l as a function of 9 and r. We see that a hole was formed, and developed 
into a dark soliton through the action of the gain-loss mechanisms. A more detailed view is 
shown in Fig. |2], which shows and the phase of a clean soliton, as a function of 9 at a fixed 
time. The signature is that the phase jumps by 7r across the soliton. The modulus does not 
precisely vanish due to small admixtures of non-soliton excitations. Fig. [3] shows the angular 
momentum per particle, which rises from to almost 1 as the condensate is being stirred, 
but settles down to ft/2 characteristic of a kink, under actions of the gain- loss mechanism. 
Fig. |] shows the total number of atoms in the ring as a function of time. The steady-state 
oscillations indicate an output train of pulses coherent with each other. In this simulation, 
the stroboscope performs the dual task of cleaning up the soliton, and acting as outcoupler 
for the mode-locked laser. In general, these tasks could be done by separate mechanisms. 
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Figure 2: Modulus and phase of the dark soliton, as functions of the angle around the ring, at 
a time when it has become fully stabilized. The phase jumps by 7r across the dip, indicating 
that the soliton is a kink. 
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Figure 3: The angular momentum went through transients while the condensate was being 
stirred, but stabilizes to a value h/2 characteristic of a kink, under actions of the gain- loss 
mechanism. 



Our numerical results are quite sensitive to the exact strength of the inter-atomic repulsive 
potential and to the gain and loss parameters. In actual experiments, however, one can adjust 
the potential, either by taking advantage of the similar couplings that exist in two-component 



Bose gases M, or by using tuning techniques involving Feshbach resonances [10], which have 



been recently demonstrated experimentally |TTf . It may also be feasible to utilize inelastic 
atom-atom collisions as a saturation mechanism, to further increase the mode-locked laser 
stability. 

We would like to acknowledge useful discussions with M. Holland. Two of the au- 
thors (A.E. and K.H.) were supported in part by a DOE cooperative agreement DE-FC02- 
94ER40818. P.D. and K.K. acknowledge ARC for the support of this work. 
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Figure 4: Mode- locking — production of a coherent train of pulses — is indicated by the 
steady-state oscillation of the total number of atoms in the condensate. 
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